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1 .0  SUMMARY 

Thib  project  is  a continuation  of  a development 
program  for  flexural  mode  transducers  and  hydroacoustic 
technology  at  Hydroacoustics  Inc.  Mathematical  models 
for  general  flexural  mode  transducers  and  hydroacoustic 
power  oscillators  developed  in  the  preceding  projects 
have  been  used  to  generate  design  trade-off  data  for 
bender  hydrophones  and  a class  of  sound  sources  which 
operate  on  the  ambient  sea  pressure  head. 

The  hydrophone  trade-off  analysis  shows  that 
bilaminar  hydrophones  are  generally  preferred  over  the 
trilaminar  configuration,  the  edge  support  for  the 
bender  can  be  designed  such  that  it  has  little  effect 
on  the  hydrophone  performance  when  operating  air  backed 
but  has  a significant  effect  on  liquid-filled,  pressure 
compensated  hydrophone  performance,  and  thicker  discs 
are  preferred  over  thinner  discs  in  the  design  of  liquid 
filled,  pressure  compensated  hydrophones.  Detailed 
design  curves  are  presented  which  quantify  the  trade- 
off between  sensitivity,  impedance,  bandwidth,  and 
operating  depth  for  many  combinations  of  hydrophone 
geometry. 

The  results  from  the  design  trade-off  study  of 
hydroacoustic  power  oscillators  describe  quantitativley 
how  the  input  and  output  power  and  harmonic  content 
change  when  the  design  parameters  valve  stroke,  supply 
pressure,  porting  length,  load  resistance,  and  filter 
characteristics  are  varied.  This  trade-off  data  enables 
the  designer  to  select  the  proper  combination  of  design 
parameters  to  give  the  desired  performance  of  the 
hydroacoustic  power  oscillator. 
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2 . 0 BENDER  H1LR0PH0NE  DES ION  TRADEOFF  STUDY 


Piezoelectric  elements  excited  in  their  flexural 
mode  have  been  used  as  hydrophones  for  the  last  20  or 
3C  years,  and  because  of  their  good  coupling  to  the 
medium,  interest  in  them  remains  high.  Flexural  mode 
(or  bender)  hydrophones  using  piezoelectric  discs 
bonded  together  or  to  a metalic  supporting  structure 
can  have  a low  characteristic  acoustic  impedance  and 
low  resonant  frequency  giving  good  power  transfer 
to  their  preamplifier  making  them  excellent  low  frequency 
hydrophones.  With  proper  pressure  compensation  this 
type  of  hydrophone  can  be  used  at  great  depth  with  no 
particular  sacrifice  in  sensitivity.  This  section  of 
the  report  presents  basic  design  and  tradeoff  data  for 
bender  hydrophones  and  shows  how  these  relationships  are 

. 

affected  by  various  edge  supports  and  pressure  compen- 
sation. Table  2.1  is  a glossary  of  symbols  and 
nomenclature  used  in  this  tradeoff  st\dy. 
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Figure  2.1  shows  the  basic  geometries  considered 
in  this  report.  The  upper  portion  of  this  figure 
shows  the  symmetrical  trilaminar  hydrophone  'where  the 
total  ceramic  thickness  t is  split  equally  on  either'"" 
side  of  the  backing  material  with  a thickness  t^.  If 
no  backing  material  is  used,  this  geometry  reduces 
to  the  simple  symmetrical  bilaminar  or  bimorph  hydrophone. 

The  asymmetrical  bilaminar  hydrophone  shewn  in  the  lower 
portion  of  this  figure  has  all  of  the  ceramic  material 
on  one  side  of  the  backing.  When  this  geometry  is 
used  with  the  water  on  the  ceramic  side,  the  hydrostatic 
pressure  will  provide  a compressive  bias  stress  in  the  |.. 

ceramic.  As  shown  later  in  this  section,  this  lets 
the  asymmetrical  bilaminar  hydrophone  be  used  to  greater 
depths  than  the  trilaminar  hydrophone. 


TABLE  2.1 

SYMBOLS  USED  IN  SECTION  2.0 

disk  radius 
disk  area  (7T a2) 
disk  total  thickness 
metal  backing  thickness 

total  ceramic  thickness 

reference  plane  coordinate 

density  of  ceramic 
density  of  metal  backing 

modulus  of  elasticity  for  the  metal  backing 
1 /S^  ceramic  material 

elastic  compliance  coefficient  at  constant 
charge  density 

elastic  compliance  coefficient  at  constant 
charge  density 

elastic  compliance  coefficient  at  constant 
charge  density 
piezoelectric  constant 

i 

dielectric  constant  at  constant  stress  (free) 
planar  coupling  coefficient 

(S?2/3^)  Poison's  ratio  of  ceramic  at  constant 
charge  density 
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TABLE  2.1  (Continued) 


(S^2/S^i)  Poison's  ratio  of  ceramic  at  constant 
voltage  (short  circuit) 


Poison's  ratio  of  metal  backing 


(Me)h  - hydrophone  receiving  sensitivity 


low  frequency  capacitance 


short  circuited,  unloaded  resonant  frequency 


Maximum  depth  based  on  maximum  radial  stress  in 


the  ceramic  material 


hinge  thickness 


hinge  outside  radius 


acoustical  compliance  of  the  disc  and  hinge 
air  backed  and  open  circuited 


acoustical  compliance  of  pressure  compensation 


cavity 


maximum  radial  tensile  stress  in  the  hinge 


maximum  radial  stress  in  the  ceramic  of  the 


bender  (compressive  stress  for  bilaminar 
geometry  or  tensile  stress  for  trilaminar 
geometry) 


ClviU  (if  PriMrt) 


Since  the  dipole  sensitivity  of  a single  flexing 
disc  is  low,  the  single-ended  hydrophones  considered 
here  are  enclosed  in  a massive  housing  which  acoustically 
terminates  one  side  in  a closed  cavity  with  stiff  walls. 
This  cavity  can  he  at  atmospheric  pressure  or  can  he 
coupled  through  a high  acoustic  impedance  to  the  ambient 
pressure  providing  pressure  compensation.  Two  flexing 
discs  can  he  used  in  one  hydrophone  to  provide  increased 
sensitivity  or  to  measure  velocity  rather  than  pressure. 
Many  of  the  results  described  here  can  he  applied  to 
both  single  and  double-ended  hydrophones. 

Some  of  the  hydrophones  discussed  in  this  report 
use  the  simple  hinge  edge  support  shown  in  Figure  2.1, 
and  other  data  assumes  a simply  supported  edge.  In 
the  latter  case  the  bender  is  assumed  to  be  simply 
supported  at  radius  a and  no  hinge  is  used.  When  air 
backed  benders  are  discussed,  the  pressure  compensation 
cavity  shown  here  is  assumed  to  have  zero  stiffness, 

. and  the  average  pressure  in  the  cavity  is  assumed  to  be 

zero. 


2.1  DESCRIPTION  OF  THE  PROBLEM 

Six  nondimens ional  design  parameters  are  defined 
from  Figure  2.1  to  describe  the  geometry  of  the  tender 
hydrophone.  Three  are  used  to  describe  the  disc:  the 
ratio  of  total  disc  thickness  to  disc  radius  (h/a),  the 
ratio  of  ceramic  thickness  to  backing  thickness  (t^/t.^), 
and  the  ratio  of  backing  thickness  to  total  thickness 
(tb/h).  These  last  two  design  parameters  are  directly 
related. 


<w 


1-(tb/h) 
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Two  parameters  describe  the  hinge:  the  ratio  of  hinge 
thickness  to  backing  thickness  (th/tb),  and  the  ratio 
of  hinge  radius  to  disc  radius  (r^/a).  The  final 
design  parameter  is  the  ratio  of  the  bender  acoustical 
compliance  including  both  the  disc  and  hinge)  to  the 
acoustical  compliance  of  the  pressure  compensation 
cavity,  (C^/C^,).  The  tradeoff  data  presented  will 
show  how  the  bender  hydrophone  performance  is  changed 
as  these  nondimensional  design  parameters  are  varied. 

Practical  design  considerations  place  general 
limits  on  these  parameters.  The  disc  thickness  para- 
meter (h/a)  can  generally  be  limited  to  a range  from 

0. C3  to  0.3.  Hydrophones  designed  with  (h/a)  below 
this  range  will  be  very  thin  and  have  high  sensitivity 
and  capacitance,  with  low  resonant  frequency,  very  shallow 
operating  depth  air  backed,  and  will  be  very  difficult 

to  pressure  compensate.  With  (h/a)  above  this  range 
the  hydrophone  will  have  the  opposite  characteristics 
from  above  and  have  its  performance  further  compromised 
by  the  effects  of  shear  displacement.  Physically  the 
parameter  (tfc/h)  is  limited  to  the  range  from  0.0  to 
1.0,  i.e.,  no  backing  to  no  ceramic.  For  the  trilaminar 
geometry  the  no  backing  case,  (t^/h)  = 0 or  (tc/t.)  = oo  , 
is  of  interest  since  this  gives  a symmetrical  bilaminar 
or  bimorph  hydrophone.  Cases  where  (tb/h)  approaches 

1 .  ^ are  of  little  interest  because  the  hydrophone 
would  have  very  little  ceramic  material.  The  range  of 
general  interest  for  the  trilaminar  geometry  is  from 
0.0  to  about  0.75  for  (tfc/h)  or  oo  to  0.33  for  (tQ/tb). 

For  the  bilaminar  geometry  these  ranges  are  from  0.20 
to  0.75  for  (tb/h)  or  4.00  to  .33  for  (tc/tb).  limits 
on  the  hinge  parameters  can  similarily  be  defined. 
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The  range  of  interest  for  the  pressure  compensation 
parameter  (CDH/CT)  is  from  0.1  to  1.0,  and  is  set  as 
much  by  performance  as  by  physical  considerations. 

Four  basic  performance  parameters  are  defined  for 
the  bender  hydrophones  considered  here:  the  low 

frequency,  open-circuited  receiving  sensitivity, 

(Me)h,  (designated  as  M for  this  report);  the  low 
frequency  electrical  capacitance,  C^q^;  short- 

circuited,  unloaded  resonant  frequency,  fR,  (designated 
here  as  f );  and  the  maximum  operating  depth,  D, 

(assuming  a safety  factor  of  about  two  for  the  ceramic 
material).  The  usable  bandwidth  of  the  hydrophone 
will  be  some  fraction  of  the  unloaded  resonant  frequency 
which  should  be  chosen  by  the  designer.  The  change 
in  sensitivity  will  be  about  3 dB  between  0.5  and  0.6 
f and  will  be  1 dB  at  about  0#4  fp* 

A nondimensional  performance  parameter  (T^/T)  is 
defined  which  relates  the  maximum  tensile  stress  in 
the  hinge  to  the  maximum  compressive  stress  in  the 
ceramic  of  a bilaminar  disc  or  the  maximum  tensile 
stress  in  the  ceramic  or  a trilaminar  disc.  With  a 
good  grade  of  steel  as  the  backing  and  hinge  material 
acceptable  values  of  this  parameter  can  be  as  high 
as  20.  for  the  trilaminar  geometry  and  2.0  for  the 
bilaminar  geometry.  When  high  strength  steel  is  used, 
these  numbers  can  be  increased  by  50#  to  30.  and  3.0. 

2.2  IDEALIZED  SIMPLY  SUPPORTED  AIR-BACKED  BENDER 

HYDROPHONES 

A useful  set  of  basic  design  curves  can  be  generated 
assuming  the  bender  hydrophone  has  an  ideal  simply  sup- 
ported edge.  These  curves  are  presented  in  parametric 


form  using  the  two  design  parameters  for  the  disc: 


(h/a)  and  (to/tb)  or  (t^/h).  The  effects  of  adding 
physically  realizable  edge  support  and  pressure  com- 
pensation are  then  expressed  as  changes  from  these 
idealized  design  curves. 


The  parameter  (h/a)  is  very  useful  for  the  design 
bender  hydrophones.  Using  it,  a family  of  hydrophones 
can  be  generated,  each  having  the  same  sensitivity, 
resonant  frequency,  or  impedance,  etc.  Thus  a large 
thick  bender  hydrophone  can  have  the  same  sensitivity 
as  a small  thin  hydrophone,  but  the  other  performance 
parameters  of  these  two  hydrophone  will  differ  greatly. 
This  gives  the  designer  a great  deal  of  flexibility  in 
specifying  his  design  to  meet  particular  performance 
goals . 


2.2.1  Basic  Design  Curves 


Figures  2.2  through  2.5  present  basic  design 
curves  which  were  derived  using  the  mathematical  model 
for  flexural  mode  transducers  developed  at  Hydroacoustics 
Inc.  under  contracts  N001 40-71  — C— 0081  and  N00024-72-C-1 295 
This  mathematical  model  predicts  the  static  and  dynamic 
stress  and  performance  of  a general  flexing  disc  with 
radial  symmetry,  axisymmetric  load,  and  sinusoidal 
excitation.  This  model  allows  for  changes  in  material 
thickness  and  properties  as  a function  of  radius  and 
does  not  rely  on  a neutral  stress  plane  concept.  Very 
general  edge  conditions  for  a flexural  mode  trans- 
ducer can  be  readily  incorporated.  This  model  also 
includes  an  approximation  for  the  shear  deflection 
which  becomes  important  for  thicker  discs  (thickness 
approaching  one  quarter  the  radius). 
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Figure  2.3  Low  frequency  capacitance  per  unit  radius 
for  bilarainar  and  trilaminar  bender  hydro- 
phones with  PZT5A  ceramic  and  steel  back- 
ing. 
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Figure  2.4  Open  circuit,  unloaded  resonant  frequency 
times  radius  for  bilaminar  and  trilaminar 
bender  hydrophones  with  PZT5A  ceramic  and 
steel  backing. 
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Firure  2.5  Estimated  maximum  depth ‘for 
trilaminar  bender  hydrophone 
ceramic  and  steel  backing. 
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The  data  presented  in  Figures  2.2  through  2.5 
are  for  a steel  backed  hydrophone  with  PZT-5A  ceramic. 

The  material  properties  used  are  shown  in  Table  2.2. 

The  equations  in  Section  2.2.3  show  how  to  modify  these 
curves  for  changes  in  ceramic  or  backing  material. 

The  performance  parameters  shown  here  are  the  same  as 
those  defined  in  Section  2.0. 

Figure  2.2  shows  that  the  receiving  sensitivity* 
per  unit  radius  is  inversely  proportional  to  (h/a). 

This  curve  also  shows  that  the  symmetrical  bilaminar 
and  asymmetrical  bilaminar  with  (t^/h)  = .4  or  (t^/tfc)  = 
1.5  have  the  highest  sensitivity.  Changing  (tb/h)  reduces 
h/a  by  as  much  as  4 dB  for  the  range  of  parameters 
shown  here.  Changing  from  steel  to  brass  or  aluminum 
backing  material  changes  M/a  by  less  than  2 dB. 

Figure  ?.3  shows  that  the  low  frequency  capacitance 
per  unit  radius  for  the  trilaminar  hydrophones*  is  about 
4 times  higher  than  that  for  the  bilaminar  configurat.'-stn 
and  that  it  is  inversely  proportional  to  (h/a). 

Changing  parameters  (tb/h)  or  (tc/tb)  over  the  ranges 
considered  here  changes  CLQW/a  by  about  6 dB.  Switching 
from  steel  to  either  brass  or  aluminum  has  a negligible 
effect  on  CLQW/a. 


*The  data  in  Figures  2.2  and  2.3  for  the  trilaminar 
hydrophones  assume  that  the  ceramic  discs  are  wired 
in  parallel.  If  the  discs  were  wired  in  series,  the 
receiving  sensitivity  would  increase  by  6 dB  and  the 
low  frequency  capacitance  would  decrease  (increasing 
the  impedance)  by  a factor  of  4. 
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Figure  2.4  shows  that  the  open  circuit,  unloaded 
resonant  frequency  times  radius  is  proportional  to 
the  parameter  (h/a)  and  does  not  significantly  change 
as  (t^/h)  or  (tc/t  ) is  changed.  The  effect  of  shear 
in  the  thicker  discs  (h/a  equal  to  .3)  reduced  f^#a 
by  5 to  15$  depending  on  the  configuration  and  value 
of  (t^/h)  or  (tc/t^).  Changing  from  steel  to  brass 
can  reduce  f^*a  by  as  much  as  15$  while  changing  from 
3teel  to  aluminum  can  change  f^*a  by  as  much  as  +10$ 
depending  upon  the  configurations. 

Figure  2.5  shows  that  changing  the  parameters 
(t,D/h)  or  (t,/tb)  has  a small  effect  on  the  maximum 
depth:  5$  for  the  range  shown  for  bilaminar  benders 
and  25$  for  trilaminar  benders.  The  maximum  depth 
shown  here  is  calculated  from  the  maximum  radial 
stress  produced  in  the  ceramic  by  a uniform  hydrostatic 
pressure  on  one  side  of  the  bender.  A safety  factor 
of  about  two  has  been  assumed  in  specifying  the  maximum 
depth. > Because  the  asymmetrical  bilaminar  bender  has 
ceramic  on  only  one  surface,  it  can  be  assembled  so 
that  a hydrostatic  load  produces  only  compressive 
stress  in  the  ceramic.  This  gives  a factor  of  ten 
increase  in  operating  depth  over  the  trilaminar  config- 
uration. This  is  reflected  in  Figure  2.5  by  two 
vertical  scales;  one  for  bilaminar  and  one  for  tri- 
laminar benders. 

2.2.2  Combined  Performance  Parameters 


It  is  possible  to  define  many  combinations  of 
the  four  performance  parameters  shown  in  Figures  2.2 
through  2 . 5 which  are  independent  of  the  parameter 
(h/a)  and  depend  only  on  the  parameter  (t^/'h)  or 
(t  /t.).  Three  such  combinations  are  described  below. 
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C1  - M * fR 


C3  "=  I#i  * CLOV  B 1 a2 

If  one  performance  parameter  is  held  constant, 
these  combinations  of  parameters  may  be  used  to  show 
the  trade-off  relationships  for  the  remaining 
parameters  as  shown  in  Table  2.3. 

Table  2.4  gives  the  values  of  the  three  combination 
parameters  for  various  values  of  (tfc/h)  and  (t0/tfe) 
for  both  bilaminar  and  trilaminar  benders.  All  three 
of  these  parameters  go  to  zero  as  (t^/h)  approaches 
zero  for  the  bilaminar  benders  and  go  to  zero  as  (t^/h) 
approaches  one  for  both  bilaminar  and  trilaminar 
benders . 


The  parameter  G^for  asymmetrical  assemblies  is 
up  to  33/>  higher  than  that  for  trilaminar  benders 
with  the  same  proportion  of  metal  backing  (t^/h). 

An  optimum  value  of  (t^/h)  exists  for  asymmetrical 
benders  whereas  this  product  decreases  monotonically 
for  increasing  (t^/h)  of  trilaminar  assemblies. 

Bilaminar  assemblies  have  values  of  C2  which 
are  approximately  4 times  that  of  corresponding 
trilaminar  assemblies.  Although  M values  are  almost 
equal  for  comparable  bilaminar  and  trilaminar  assemblies, 
the  values  of  are  not.  Parameter  O2  bas  a 

maximum  value  near  (t^/h)  equal  to  .2,  (t  /tfa)  equal 
to  4.,  for  bilaminar  benders.  No  such  peak  exists 
for  the  trilaminar  configuration. 
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Table  2.3 

Design  Tradeoff  Parameters  For 
One  Performance  Parameter  Fixed. 


Design  Criteria 


Performance 

Parameter 


Fixed  Depth 


Fixed  Sensitivity  Fixed  Fixed 

Capacitance  Bandwidth 


Table  2 .4 


Combined  Parameters  for  Bilaminar  and  Trilaminar 
Bender  Hydrophones  with  PZT-5A  Ceramic 
and  Steel  Backing 


Trilaminar 


,(1) 


tc/tb 

tb/h 

C1  C2 

C3 

oc 

0.0 

31.5  0.259 

6.17 

2.33 

0.3 

27.5  0.151 

8.68 

0.67 

0.6 

18.8  0.057 

10.33 

0.00 

1 .0 

0.0  0.0 

0.0 

Bilaminar 

tc//tb 

Vh 

°1  C2 

C3 

00 

0.0 

0.0  0.0 

0.0 

4.00 

0.2 

32.7  0.624 

20.0 

1 . 30 

0.4 

39.2  0.585 

32.9 

0.67 

0.6 

33.8  0.356 

42.6 

0.00 

1.00 

0.0  0.0 

0.0 

* 

se 

ii 

\ — 
O 

r 

' R 

1C"10*  KHz  * 

V/  MPa 

O 

II 

CL0W 

1 0" 1 3*  V/  (>aPa 

* pfd) 

C.J.  - M * 

CL0W  * 1 * * / a 

10“4  * V * pfd  * ft/ (MPa 

(1)  The  data  for  the  trilaminar  hydrophones  assume  that 

the  ceramic  discs  are  wired  in  parallel.  If  the 

discs  are  wired  in  series,  the  sensitivity  would 
increase  by  6 dB  and  the  low  frequency  capacitance 

would  decrease  by  a factor  of  4. 


The  parameter  shows  that  in  general  sensitivity 
must  he  traded  off  for  depth  unless  capacitance 
can  he  reduced.  Parameter  is  greater  for  bilaminar 
assemblies  because  of  their  greater  depth  capability 
when  assembled  to  be  compressively  stressed  by 
hydrostatic  pressure.  This  parameter  ha3  a maximum 
value  above  (tb/h)  equal  to  .6,  (t  /tfe)  equal  to  .67, 
for  both  the  bilaminar  and  trilaminar  benders. 


2.2.3  Equations  to  Account  for  Material  Changes 

This  section  presents  equations  which  can  be 
used  to  approximate  the  change  in  performance  para- 
meters when  either  the  ceramic  or  backing  material 
properties  are  changed.  These  equations  are  simplified 
where  possible  and  use  the  nomenclature  defined  in 
Table  2.1.  One  simplification  is  made  with  respect  to 
the  Poisson  ratio  for  the  backing  and  ceramic.  For 
most  metals  of  interest  is  approximately  0.3, 

and  for  most  ceramics  U3ed  in  hydrophones  o'®  is 
approximately  0.5.  Parameters  which  are  weak  functions 
of  the  Poisson  ratio  use  these  approximations. 

2 . 2 . 3 • 1 Symmetrical  Trilaminar  Equations 

The  equations  for  the  performance  parameters 
of  a symmetrical  trilaminar  bender  hydrophone  are 
directly  related  to  the  assumed  normalized  deflection 
curve  of  the  bender.  These  equations  are  developed 
from  the  closed  form  solution  of  a symmetrical 
trilaminar  bender.  Four  coefficients  can  be  defined 
which  are  independent  of  this  deflection  curve;  five 
other  coefficients  must  be  evaluated  as  a function 
of  the  bender  shape. 
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(h  4 tfe) 

4 tvh  4 t; 


\>H  + 


For  voltage  parameters  M,  and  C^q^. , tne  assumed 
normalized  deflection  curve  is 


h = 1 


- M 


which  is  the  curve  for  an  applied  edge  moment. 

Because  of  its  simplicity,  this  curve  is  also  used 
to  approximate  the  maximum  radial  stress  which  in 
turn  is  used  to  estimate  the  maximum  depth.  The 
following  five  coefficients  were  developed  using  this 
deflection  curve.  / \ i /? 


K = 27Tkp(h  4 tb) 


20  - 

0V  /a  r \ 

s1 1 (i  - <r  ) 


b _ 1 4-  cr  ^ 


*B  i 4<rI)  = .75 

/\  = r} 
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The  deflection  equation  used  to  calculate 
; short  circuit,  unloaded  resonant  frequency  s 

Y = 1 - 1.4417(|)  2 + -”04  (|)  3 1-  -HI2  ( 
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and  the  following  coefficients  are  defined. 


= .353  + .262  <rl  = .432 


Ar  = .353  + .262  <rD  = .484 

Ae  = AE  - .231  k£  (1  -a-A/C,  = .484  - .116  k*  < 


16.8 

T" 

a“ 


1 .33  A 


E,v3 


W-*Z) 


+ .48  Ebt 


2. 2. 3. 2 Asymmetrical  Bilaminar  Equation 

A simple  neutral  stress  plane  development  of 
the  model  for  asymmetrical  "bilaminar  discs  does  not 
yield  satisfactory  prediction  accuracy.  The  fol- 
lowing equations  have  been  empirically  developed 
from  the  computer  model  predictions  to  predict  the 
approximate  performance  for  the  t^/h  ratio  which 
gives  the  maximum  receiving  sensitivity. 

A coordinate  of  the  mechanical  neutral  stress 
plane,  t , has  proven  to  be  a useful  parameter  in 
locatingnthe  near  optimum  ratio  of  ceramic  to  metal 
thickness  and  is  defined  below. 
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A positive  value  of  tn  means  that  this  plane  is  in 
the  hacking  material;  a negative  value  places  it  in 
the  ceramic.  The  optimum  value  of  t^/h  is  slightly 
greater  than  or  equal  to  that  obtained  when  tn  - 0.0 


j?  - 

(?/Ec  - 
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The  approximate  equations  for  the  performance 
parameters  for  this  value  of  t^/h  are  given  below 

2 g^ 

kp  ° 2 ~ S^(1  VV£33 


cb  = A(1  -kp)£33/tc 


I 


CL0W  = 1.47  C 
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2.3  AIR-BACKED  HYDROPHONE  WITH  HINGE  EDGE  SUPPORT 


In  the  previous  section  computer  generated 
design  trade  off  curves  for  bender  hydrophones  were 
presented  as  functions  of  two  non-dimensional  para- 
meters: (h/a),  total  bender  thickness  to  radius 

ratio,  and  (*c/t^)>  ceramic  thickness  to  backing 
thickness  ratio.  A third  parameter  (t^/h)  which  is 
directly  related  to  (t  /tb)  was  also  shown.  This 
data  assumed  the  simplest  bender  hydrophone  geometry: 
a bender  with  uniform  thickness  for  both  the  backing 
and  ceramic  material  and  an  ideal  simply  supported 
edge.  The  performance  parameters  used  for  these 
curves  are  open  circuit  receiving  sensitivity,  M; 
low  frequency  electrical  capacitance,  CL0W;  short- 
circuited,  unloaded  resonant  frequency,  fpj  and 
maximum  depth,  D.  The  effect  on  these  parameters 
of  adding  a simple  hinge  to  the  bender  hydrophone  is 
now  examined. 

2.3.1  Simple  Uniform  Hinge 

A sketch  of  the  composite  bender  and  hinge  to  be 
considered  here  was  shown  in  Figure  1 in  Section  2.0. 

Two  non-dimensional  parameters  are  used  to  describe 
the  hinge:  hinge  thickness  to  backing  thick- 

ness ratio;  and  (r,n/a),  hinge  radius  to  ceramic  radius 
ratio.  One  additional  performance  parameter  has  been 
introduced  to  describe  the  stress  in  the  hinge:  (r^/a  ), 
maximum  radial  hinge  stress  to  maximum  radial  ceramic 
stress  ratio.  Most  backing  materials  can  support 
twice  the  compressive  stress  and  twenty  times  the 
tensile  stress  possible  for  the  ceramics  used. 


i 
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The  maximum  radial  bending  stress  in  the  hinge  is 
strongly  dependent  on  the  hinge  thickness.  For  a 
uniform  static  pressure  this  stress  is  roughly  pro- 
portional to  for  between  0.25  and  1.0 

(for  trilaminar  bender  parameters  (h/a)  = 0.03  and 
(t  /t,  ) = 2.33).  This  stress  would  be  proportional 

to  t^  "'for  an  ideal  clamped-c lamped  ring. 

When  a simple  hinge  like  that  shown  in  Figure  2.1 
is  used  instead  of  an  ideal  simple  support,  the 
parameters  Cp^,  fR,  and  D change  as  shown  in 
Tables  2.5  and  2.6.  These  results  apply  to  any  value 
of  the  relative  thickness  parameter  (h/a)  except  for 
fp  for  which  the  change  factor  is  valid  for  thin  hinges 
but  is  only  an  approximation  for  thicker  hinges  because 
of  the  effects  of  shear. 

These  tables  are  used  with  the  four  basic  design 
trade  off  curves  presented  in  the  previous  section. 

To  assess  the  effect  of  adding  a simple  hinge  to  a 
bender  hydrophone,  multiply  any  of  the  performance 
parameters  K,  fR,  or  D by  the  appropriate  factor 

shown  in  this  table.  The  last  entry  in  this  table 
shews  the  ratio  of  hinge  to  ceramic  stress.  As 
mentioned  earlier,  values  of  this  parameter  greater 
than  2.0  for  bilaminar  and  20.0  for  trilaminar  hydro- 
phones may  give  a design  in  which  the  hinge  stresses 
limit  the  operating  depth.  All  of  the  bilaminar  hydro- 
phone designs  shown  here  have  a stress  ratio  greater 
than  2.0  although  some  of  them  are  close  enough  to 
be  considered.  (Materials  with  an  ultimate  tensile 
strength  greater  that  200K  psi  would  be  required  for 
the  backing  material.) 


Changes  in  Performance  Parameters  When 
A Simple  Hinge  Replaces  the  Ideal  Support  for 
Air  Backed,  Eilaminar,  Bender  Hydrophones 


rh/a  = 1.10 


rh/a  =1.20 


(W 

1.00 

.75 

.50 

1 .00 

.50 

(tc/tb) 

H 

4.00 

1.50 

1.53 

1.52 

1.78 

1 .74 

1.50 

.885 

1 .09 

1.24 

1 .-23 

1 .50 

.67 

.425 

.695 

1 .03 

.750 

1.35 

.55 

.251 

.470 

.848 

n 

4.00 

. 9^2 

.972 

. 979 

.973 

.984 

"LOW 

1.50 

.872 

.926 

.970 

.918 

.982 

.67 

.807 

.857 

.928 

.843 

.956 

.55 

.348 

» 

.874 

.924 

4.00 

.92$  ' 

.888 

.788 

.759 

.421 

iR 

1.50 

1.04 

.958 

.887 

.914 

. 727 

.67 

1 .22 

1 .08 

.931 

1 .08 

.797 

.55 

1 .31 

1.16 

.968 

D 

4.00 

1.07  ‘ 

1 .0$““ 

1 .02 

.942 

.900 

1 .50 

1.19 

1 .06 

.986 

.988 

. 863 

.67 

1 .56 

1 .27 

1.03 

1 .22 

.864 

.55 

1 .91 

1.52 

1.12 

T-u/T 

4.00 

2.37 

3.54  ‘ 

7.71 

5.32 

20.3 

h' 

1.50 

2.19 

2.11 

2.51 

2.22 

5.10 

. 67 

2.70 

2.77 

2.62 

2.33 

3.05 

.53 

3.03 

3.44 

3.36 

- 

— 

Table  2.6 


Changes  in  Performance  Parameters  When 


A Simple  Hinge 

Replaces 

the  Ideal 

Support  for 

Air  Backed, 

Trilaminar 

■,  Bender 

Hydrophones 

Va  = 

1.10 

rh/a  =1.20 

<W 

1.0 

0.5 

1.0  0.5 

M 

2.33 

.851 

1.15 

1.19  1.40 

.67 

- 

- 

1 .22 

CL0W 

2.33 

.67 

.879 

.978 

.924  .988 

.957 

fR 

2.33 

.994 

.833 

.822  .655 

.67 

- 

- 

.797 

P 

2.33 

1 .09 

.924 

.906  .817 

.67 

- 

- 

.882 

These  results  show  that  thin  hinges  with  larger 
values  of  (r^/a)  would  be  preferred  if  the  stress 
criterion  can  be  met.  In  general  these  hinges  would 
slightly  increase  the  hydrophone  sensitivity  over  an 
ideal  simple  edge  support  and  have  a small  effect  on 
the  other  parameters. 


The  low  frequency  capacitance,  is  always 

reduced  by  the  addition  of  a simple  hinge.  The 
resonance  frequency  fp  was  also  reduced  for  most  of 
the  hinge  configurations  tested.  The  results  also 
show  that  in  general  the  depth  limit,  I),  of  the  air 
backed  bender  hydrophone  is  reduced  with  a clamped 
hinge  design.  Increases  in  the  depth  capability  of 
the  hydrophone  gained  by  the  addition  of  a simple 
hinge  are  always  accompanied  by  a reduction  in  the 
sensitivity.  Table  6 is  a much  abbreviated  version 
of  Table  5 but  shows  the  same  trends  for  the  trilaminar 
bender  hydrophone  as  seen  for  the  bilaminar  geometry. 


2.3.2  Tapered  Web  Hinge 


One  practical  means  of  approximating  a simple 
uniform  hinge  is  shown  in  Figure  2.6  where  a flexible 
edge  support  is  achieved  by  a smoothly  tapered  web. 

The  performance  of  this  hinge  is  determined  using 
Table  2.7  and  Figures  2.2  through  2.5  in  a way  similar 
to  that  described  in  the  previous  section.  Each  change 
factor  relates  the  performance  of  a hydrophone  with 
this  web-type  hinge  to  the  comparable  hydrophone 
with  an  ideal  simple  edge  support. 
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CHANGES  IN  PERFORMANCE  PARAMETERS  WHEN  A WEB-TYPE 
HINGE  REPLACES  THE  IDEAL  SIMPLE  SUPPORT  FOR  AIR- 
BACKED,  BI LAMINAR,  BENDER  HYDROPHONES 


rh/a  = 

1.10 

Vs 

i 
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.25 

.125 

.0625 
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Figure  2.6  Geometry  and  Dimensions  of  Tapered  Web  Hinge 


For  this  study  the  ratio  of  minimum  hinge 
thickness  to  back  material  thickness  (t^/t.  ) was 
varied  from  0.5  to  0.0625.  The  performance  with  this 
web  hinge  design  is  about  the  same  as  that  found  for 
a simple  uniform  hinge.  The  data  shows  a slight 
reduction  in  the  stress  ratio  (T^/T),  but  it  is 
possible  that  the  true  maximum  stress  wa§  not  found 
in  the  web  hinge.  The  region  of  highest  tensile 
stress  occurs  at  a radius  midway  between  the  radius 
of  minimum  and  maximum  thickness  in  the  clamped  edge 
side  of  the  hinge.  The  stress  data  shows  that  webs 
as  thin  as  0.125  t^  probably  are  not  useful. 


2.4  PRESSURE-COMPENSATED  HYDROPHONE  WITH  SIMPLE  HINGE 

EDGE  SUPPORT 

The  hydrophones  described  up  to  this  point  have 
been  limited  to  relatively  shallow  water  'down  to 
several  thousand  feet).  Many  applications  require 
operation  in  water  with  depths  greatly  exceeding  2000 
feet.  The  hydrostatic  pressure  at  these  depths  will 
fracture  the  disc  unless  some  form  of  pressure  compen- 
sation is  provided.  Since  <the  performance  of  a bender 
hydrophone  can  be.  affected  by  the  impedance  of  the 
terminating  cavity,  an  additional  design  compromise 
must  be  made  between  the  hydrophone  size  and  its 
acoustic  performance. 

* 

The  most  likely  (and  at  firs$  appearances,  best) 
choice  would  be  to  fill  the  terminating  cavity  with 
a gas  which  is  maintained  at  the  skme  pressure  as 
the  ambient  hydrostatic  pressure.  .Unfortunately  such 
pressure  compensation  systems  have  proven  to  be 
unreliable  and  quite  bulky  making  it  necessary  'jo 
go  to  a liquid  filled  pressure  compensation  cavity. 

A liquid  with  high  compressibility  (like  silicone  oil) 
can  be  used  in  the  cavity,  and  a pressure  equalizing 
system  with  high  acoustic  impedance* (a • long,* small 
diameter  hole)  can  be  used  to  compensate  for  the 
hydrostatic  pressure.  This  section  ’treats  such  a 
pressure  compensation  system  for  bender  hydrophones. 

The  liquid  filled  cavity  will  change  the  acoustic 
performance  of  the  hydrophone  by  adding  stiffness  which 
increases  the  resonant  frequency,  decreases  the  receiv- 
ing sensitivity,  and  decreases  the  low  frequency 
capacitance.  The  change  in  low  frequency  capacitance 
is  small  and  can  be  neglected  entirely  for  suitable 
cavity  coip)liance. 
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The  acoustic  impedance  of  the  pressure  equalizing 
passage  will  introduce  a low  frequency  roll  off  to 
the  hydrophone  frequency  response.  The  compliance 
of  the  terminating  cavity,  C^,,  is  given  by  = V/0  , 
where  V is  the  volume  of  the  liquid  filled  cavity 
and  0 is  the  bulk  modulus  of  the  liquid.  The  roll-off 
frequency  is  given  by 


1 

fLOW  = 2TWLP°T 

where  Lp  is  the  passage  inertance  which  compensates 
for  the  hydrostatic  pressure.  This  roll-off  can  be 
made  to  start  at  an  arbitrarily  low  frequency  by 
adjusting  the  passage  inertance  (changing  the  area 
and  length). 

The  effect  of  the  liquid  filled  cavity  on  the 
hydrophone  performance  can  be  described  using  the 
ratio  of  the  compliance  of  the  air-backed,  composite 
disc-hinge,  Cru,  to  the  compliance  of  the  terminating 
cavity,  Cp.  The  liquid  backed,  pressure  compensated 
hydrophones  of  most  interest  have  less  than  3*Cm 
or  the  terminating  cavity  is  less  than  three  times  as 
stiff  as  the  composite  disc-hinge.  The  sensitivity 
and  resonant  frequency  of  these  pressure  compensated 
hydrophones,  ;!P~,  and  Fpp-,,  are  related  to  the  air 
backed  sensitivity  and  resonant  frequency  by  the  follow- 
ing relationships. 


M ( 1 + °DH/CT  ) 
f R ^ 1 + ^H^T  ^ 


The  values  for  M and  f ^ are  obtained  using  the  pro- 
cedures outlined  in  Sections  2.2  and  2.3.  As  can 
be  seen  in  these  two  equations,  the  sensitivity- 
bandwidth  product  for  a pressure  compensated  hydro- 
phone is  always  less  than  the  product  for  the  same 
hydrophone  air-backed.  To  minimize  the  effect  of 
liquid  tacking,  large  liquid  compliances  are  desirable. 
The  compliance  Cn  is  proportional  to  the  cavity 
volume  and  is,  therefore,  limited  by  the  physical 
size  the  designer  is  willing  to  use.  The  cavity 
volume  is  also  limited  by  the  highest  frequency  to 
be  measured  since  the  cavity  dimensions  should  be 
a 3mall  fraction  of  a wavelength. 

A satisfactory  liquid  compensated  hydrophone 
is  not  necessarily  obtained  by  simply  adding  a 
liquid  cavity  to  a good  air-backed  hydrophone.  A 
sensitive  air-backed  bender  hydrophone  will  generally 
have  a small  (h/a)  ratio  and  be  limited  tc  quite 
shallow  operation.  When  this  hydrophone  is  liquid 
backed  to  provide  pressure  compensation,  the  results 
are  discouraging  at  best  because  a very  large  pressure 
compensation  cavity  is  required.  Since  the  basic 
compliance  of  a disc  is  inversely  proportional  to 
the  ratio  (h/a)  cubed,  thicker  discs  can  use  much  more 
reasonable  liquid  compensation  cavities.  One  disadvan- 
tage in  using  thick  discs  is  the  added  displacement 
due  to  shear  which  can  increase  the  compliance  of  the 
hinge  and  disc  by  as  much  as  40$.  The  optimum  value 
of  (h/a)  seems  to  be  near  0.1  where  the  increase  in 
sensitivity  for  a thinner  disc  is  off  set  by  the  re- 
duction in  sensitivity  due  to  liquid  backing  and  the 
shear  displacement  is  not  significant.  For  these  reasons 
thicker  discs  with  stiffer  hinges  are  more  suited  for 
pressure  compensation  than  thin  discs  or  flexible  hinges. 


Table  2.8  lists  the  acoustical  complance  for 
several  combinations  of  the  discs  and  hinges  considered 
earlier.  This  table  shows  the  importance  of  the  hinge 
when  a hydrophone  is  pressure  compensated  using  liquid 
backing.  In  Section  2.3,  it  was  shown  that  adding  a 
properly  designed  hinge  to  a disc  has  a small  effect 
on  the  hydrophone  performance  parameters.  Adding  a 
hinge  to  a disc  does  make  a significant  change  in  the 
compliance  of  the  composite  disc-hinge.  A two  to  one 
change  in  such  as  shown  in  this  table  will  require 
a two  to  one’change  in  the  volume  of  the  pressure  com- 
pensation cavity.  The  data  shown  here  are  for  bilaminar 
discs  with  a radius  of  0.01  M (1  cm)  and  a thickness 
ratio  (h/a)  of  0.3.  In  order  to  use  this  table  for 
other  values  of  a or  (h/a)  the  following  relationship 

should  be  used. 

CDH(a.  h/a)  = (100  a)3  CDH(,Cr’  0,3) 

The  values  of  CT(.01,  0.3)  are  from  the  table  and  a 
is  In  meters.  Since  the  effect  of  shear  displacement 
is  included  in  the  data  for  this  table,  the  above 
relationship  is  not  accurate  for  small  values  of  ( h/a). 
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3.0  MODELING  DIRECT  MODULATION  OF  SEA  WATER  FLOW 


The  hydraulic  valve  which  modulates  the  flow 
between  the  high  pressure  source  and  low  pressure 
sink  is  the  heart  of  the  hydroacoustxc  power  oscil- 
lator. Typically  it  has  been  analyzed  using  Bern- 
oulli's Equation*  for  the  steady  flow  of  an  inviscoid, 
incompressible  fluid. 

£3  + il  + gZ.  _ Ez  + Iz.  + SZ2 
r % - f x 

Where  P is  the  pressure,  V is  the  particle  velocity, 

Z is  the  height  above  datum,  p is  the  fluid  density, 
and  g is  the  acceleration  due  to  gravity.  Suoscripts 
1 and  2 refer  to  a section  of  fluid  on  opposite  sides 
of  the  valve.  Let  subscript  1 be  used  for  conditions 
in  the  supply  reservoir  where  V is  approximately 
equal  to  zero,  and  2 be  used  for  conditions  at  the 
exhaust  of  the  valve  which  is  at  the  same  datum 
level  as  the  supply  reservoir  (Z1  = Z2). 

Bernoulli's  equation  reduces  to 


P,  = 


P2  + 


or 


- — (P  - 

This  can  be  rewritten  as 


p2) 


Q|Q|  = 


2 2 
2A2  C2 


b? 


Where  Q is  the  volume  flow  rate,  A is  the  open 

area  of  the  valve,  and  is  the  discharge  coefficient 

(Cjj  can  be  as  high  as  1.0  for  fully  rounded  valves, 

and  C-p,  as  low  as  0.6  for  sharp  edged  valves). 

itpor  example , see  TrFIuIT  t’ower  Control"  by  J.E.  Blackburn, 
G.  Reethof , and  J.O.  Shearer,  John  Wiley  & Sons,  I960 
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The  shape  of  the  valve  opening  has  very  little  effect 
°n  Cpj  a long  narrow  slit  has  the  same  value  as  a 
circular  hole.  To  reiterate,  Equation  1 applies  to 
an  ideal  nonviscous  and  incompressible  fluid  flow 
where  steady-state  conditions  prevail.  These  assump- 
tions are  reasonably  accurate  for  most  hydraulic  valves 
used  in  fluid  power  control  systems  but  must  be 
reexamined  when  considering  the  hydroacoustic  power 
oscillator  which  operates  at  considerably  higher 
frequency  than  most  fluid  power  systems. 
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3.1  differential  equation  for  valve  flow 

Navier-Stokes  equations*  which  describe  the 
conservation  of  momentum,  constitute  the  fundamental 

basic  of  fluid  flow  theory.  If  it  is  assumed  that 
the  fluid  is  nonvisoous  or  that  the  contribution  oy 
viscous  terms  in  the  equations  is  much  smaller  that 
that  by  the  inertia  terms,  one  can  ignore  the  affects 
of  viscosity.  Then  the  equations  are  reduced  to  Euler  s 
equations  valid  for  an  ideal  fluid  and  are  given  by 

pH  + (v-  . V + P = F 

where  P is  the  fluid  density,  AT  is  the  particle  velocity 
vector,  P is  the  fluid  pressure  and  F is  the  applied 
force  vector. 

For  one  dimension  flow  this  equation  reduces  to 

0- 

and  can  be  integrated  along  a streamline  from  region Q, 
to  region  Qj  as  shown  in  Figure  3-1  . 

o-fa 


ds  +\  ds 


= P2  - p1  +m f ds  (s*n 


where  sgn  v equals  1 if  v is  positive  and  equals  -1 
if  v is  negative.  The  velocity  in  region <35  is  assumed 
to  be  zero,  and  v,  is  the  discharge  velocity  from 
the  orifice.  The  conversion  from  particle  velocity  to 
the  volume  velocity  through  the  valve,  Qv,  la  given  by 

♦For  example , see  "Mechanics  of  Deformable  Bodies" 
by  Arnold  Sommerfeld,  Academic  Press,  1950. 
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Figure  3.  1 Ceometr, 


for  Valve  Flow 


Approximate  y in  region  1 by  -10X  , by  0 in  region 

o 

2,  and  carrying  out  the  indicated  integration  gives 


SL.  dy  = 


2K 


CDXS  f 

2 J 


t 


In 


/p  V ^ 

UD  AS 


(v 


ioo  x: 
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The  square  root  terms  in  this  equation  can  be 

approximated  using  the  series  expansion. 

. 2'»  _ — 

o I 

r2  . J-r>Xf 


100X“  + 


= 10X„  1 + 


D 

800_. 


The  equation  for  the  integral  term  becomes 

dy 


aquation  for  ti 

f A 

K i » ! l\ 


where  1 + 


k ln(40/C^  ^ q _ p$v 

2k ft  210^ 


'D 


800 


has  been  approximated  by  1 . Since 


Cp  is  approximately  equal  to  0.6, 


In  (40/Cjj)  = 

2 


2.0 


and  the  final  form  of  the  differential  equations 
for  the  alternating  flow  through  the  valve  is 


£2  - K 

dt  ~ TJT 


?1 


- P, 


9 Qi 

fee 


dXS  j>Qv|  QV  I 
dt  "2K?C2X2 

dX, 


dQ  . U/VS  ... 

and  —£~r  are  negligible 


When  the  dynamic  terms , 
this  reduces  to  the  form  of  Bernoulli's  equation 
shown  at  the  beginning  of  this  section. 


3 . 2 SUPERIMPOSED  ALTERNATING  AND  STEADY  FLOW 


In  a transducer  utilizing  direct  modulation 
of  sea  water  flow  to  produce  an  acoustic  signal  in 
the  water,  the  pressures  and  flows  have  both  steady 
and  alternating  components.  When  these  pressures 
and  flows  encounter  an  orifice,  the  amount  of  pressure 
drop,  power  loss,  and  harmonic  content  all  depend  on 
the  relative  magnitude  of  the  steady  and  alternating 
components . 

3.2.1  Pressure  Drop  and  Fewer  Loss  in  an  Orifice 
When  0=QQ+01  sin  (wt)  + Q2  sin  (2  wt) 

Returning  to  the  basic  Bernoulli  equation  for 
the  flow  through  an  orifice,  one  can  find  the  pressure 
drop,  p,  and  instanteneous  power  loss,  W,  using  the 
following  equations^ 


- .V— - -W-., 
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Where?  is  the  fluid  density,  Q is  the  instantaneous 
flow  through  the  orifice,  CD  is  the  discharge 
coefficient  (CD  = .6),  and  A is  the  orifice  area. 

In  many  cases  of  hydroacoustic  devices, 
there  are  fixed  orifice  like  constrictions  which 
must  pass  both  steady  and  alternating  flow. 

Examples  are  feed  tubes  and  inertance  tubes.  In 
such  cases  the  alternating  flow  causes  an  increase 
in  the  average  pressure  drop  over  that  which 
would  occur  in  the  case  of  steady  flow  only. 

Also,  in  such  cases  the  total  power  loss  is  of 
interest  in  computing  overall  efficiency  and  power 
balance.  The  following  develop  expressions  for 
both  these  quantities.  Let  the  flow  be  given  by 

q = + Q1  sin  (wt)  + Q2  sin  (2  wt). 

i 


where  K and  M are  functions  of  and  Q2  and  Pg,p  and 

W are  the  pressure  drop  and  power  loss  respectively 
RT 

in  the  absence  of  alternating  flow. 


If  the  instantaneous  pressure  and  power  loss  are 
integrated  over  one  cycle,  the  average  pressure 
and  power  loss  are  given  by 

„ 2 


2 C 


w . vv  = M £—> 

hVE  2 C-.  A 


= " WST 


f 


If  k1 


Q1/Q0  and 


i 


i 


i 


k"  = Q2/Q1 

then  K and  M can  he  found  from  Figures  3.2  and  3.3. 

If  kn  is  equal  to  zero  and  k1  i3  large  (k,>5),  then 
K and  M can  he  approximated  hy 

K»(1  +8k,)  / 2 7T 

M^(1  + 12  k1  + 8 k1  ^/  ) /2TT. 

If  k.  is  very  large,  these  expressions  heccme 

KX4  k1  /TT 
H»4k3  /TT 

3.2.2  Harmonics  of  Pressure  with  Q = sin  (wt) 


In  the  above  section  the  pressure  factor  K is 
calculated  to  determine  the  average  pressure  difference 
for  combined  steady  and  alternating  flow  through  an 
orifice  like  constriction  when  the  average  pressure 
difference  for  steady  flow  alone  is  known.  In  this 
section  an  estimate  of  the  magnitude  and  phase  of  the 
higher  harmonics  is  derived  for  a known  flew 

Q = Qc  + Q1  sin  wt. 

The  flow  is  a result  of  a pressure  difference 


q|  Q|  where 

r 


ci  * 


* — Y 

2 C2  A2 


and 


C-j.  is  the  discharge  coefficient 


A is  the  area  of  the  orifice. 


I 
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Power  Factor,  W,  for  Supenmpo 
Alternating  and  Steady  Flow. 


Assuming  that  quasi-static  conditions  exist  for  the  flow 
field,  the  harmonic  analysis  of  P,  is  described  as 
follows.  Let 


sin  (Nuit  + i>N) 


where  is  the  magnitude  and^is  the  phase  cf  then’K 
harmonic  of  the  normalized  pressure.  The  pressure  has 
been  normalized  with  respect  to  the  pressure  in  the 
absence  of  the  alternating  flow. 


l't  can  be  shown  that  for  n>2  and  (k^=  Q<|/Qq) 

if  n is  even 

(l+ki2)  sinn01  _ [cos(n-1  )01  _ cos(n+1)0^- 
"2  n 1 1;  ~ 


An ' t 


k, 

i 

4 


sin(n-2)01  + sin(n+2)0^ 


n-2 


nh2 


if  n is  odd 


An  = it 


' - +k1  j cos  n01  + k sin(n-1)01  _ sin(n+1)0^!_ 
2 n 1 L rTi  n+1  J 


cos(n+2)01  + cos(n-2)0. 
n+5  n-2  “ 


49 


I 


A.  = 2k.  4 4 6+  kf)  2008  *>1  - ki  f_  ,fl  , 8ln  2|Z1 , 

1 f |V  — ? T (r~2V  nTS— > 


2 r 

2 = ~k1  + 4 0+k.,2)  sin  201  k fcos  0 cos  3 0-,^ 

4 ? — ’ 1 —3-/ 

+ Hi!  r^-0,  - sip  4 gi\ 

4 v * 1 4 ' 

where  01  = S^n  (k.^  and  ^=0  if  N is  Odd,  or 


= JL  if  n is  even. 
n 1 


-f  k1  - 1 , 

A1  = 2 k1 


A2  = '!i 


Ar  = 0 f or  n >2 


For  very  large  values  of  k, , 


A = 3 if  n is  odd  and 

n K n(n2-4) 

A = 0 if  n is  even. 


These  last  expressions  are  the  Fourier  series  for  k^ 
sin  wt|sin|wt  and  represent  the  limit  when  the  steady 
flow  is  zero. 


Figure  3.4  is  a plot  of  An  as  a function  of 
k1  for  n = 1 to  4.  As  can  be  seen  from  the  graph, 
the  first  harmonic  is  predominant  and  the  higher 


First  Four  Harmonics  in  the  Pressure 
Acroas  an  Orifice  when*  the  flow  is 
Q0( 1 - k,  sin  wt) . 


harmonics  keep  on  becoming  less  significant  succes- 
sively. Prom  the  above  equations,  in  the  limit  of 
large  k^ , the  magnitudes  of  the  third  and  fifth 
harmonic  are  respectively  only  about  20 % and  3 i°  of 
the  magnitude  of  the  first  harmonic  while  the  second 
and  fourth  harmonics  tend  to  dissipate  altogether. 

3.2.3  Harmonic  Content  of  Flow  Through  a Fixed  Orifice 
when  P=  sin  wt 

Since  a fixed  orifice  is  a nonlinear  element, 
it  generates  harmonics  in  the  flow  through  it  when 
alternating  pressures  are  applied  across  it.  When 
the  amplitude  of  the  oscillatory  pressure  is  large, 
the  nonlinearities  may  be  significant  and  the  higher 
harmonics  may  not  be  negligible  when  compared  to  the 
fundamental.  It  is  useful  to  know  the  relative 
magnitudes  of  the  harmonics  generated  by  this  type  of 
orifice . 

Let  the  pressure  upstream  of  the  orifice  be 

_ A 

?1  - F1  + P1  sin  wt 

and  the  pressure  downstream  be  a constant 

If  quasi-steady  state  conditions  are  assumed  to 
exist,  the  resulting  flow  is  given  by  Bernoulli's  equation. 

Q = C H F(t)  |j*  Sgn  F(t ) 

Where  F(t)  = [(^1-P2)+P1  sinUJtl 
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Then  Q = C|P,  - ?2|  * | 1+  ^sin  wt  p Sgn  (1  + j,  sin  v^t)  Sgn  (P.,- 


Where  QU  P / + Pg)  or 


^n=  £ = |l+  <*  sir.cjt  Sgn  (1+<A  sini«>t) 

Cl^-Pg]  fSgn  (?rP2) 


Where  Q is  the  flow  normalized  with  respect  to 
the  flow  caused  by  steady  pressure  difference  alone. 


A computer  program  was  developed  to  perform  a Fourier 


analysis  on  the  expression  for  Qn  and  determine  the 
magnitude  ar.d  phase  of  the  first  nine  harmonics. 


Let  Qr  = ^ Ck  cos  (kwt  +^k) 


where  is  the  magnitude  of  the  ktb  harmonic  and 


>k  is  its  phase  angle. 


* °1 


?[v, 


cos  (kwt  + 


**>] 


Qn  = C-  £ 008  Ocwt  + 


where  c£  is  magnitude  of  the  ktb  harmonic  normalized 
with  respect  to  that  of  the  first  harmonic. 


I: 


i 
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Figure  3.5  shows  the  magnitude  and  phase  angle 
of  the  first  four  terms  in  this  expression  for  various 
values  of  the  modulation  factory  . For  the  range  of 
cC  covered  here,  the  second  harmonic  is  always  in 
phase  with  the  fundamental  and  reaches  a peak  of 
.32  at  a modulation  factor  of  1.2.  The  third  harmonic 
is  less  than  .15  and  is  in  phase  with  the  fundamental 
up  to  a modulation  factor  of  approximately  2.0  at 
which  point  it  becomes  out  of  phase  with  the  fundamental. 
AsdLis  reduced  to  zero,  all  the  harmonics  approach 
zero  and  the  flow  is  given  by  the  steady  flow  component. 

It  can  also  be  seen  from  the  figure  that  as cHincreases, 
the  second  harmonic  tends  to  decrease  slowly  while  the 
third  harmonic  does  not  follow  this  pattern.  In  the  limit 
of  dt very  large,  (1  + 0(|Sin  wt  | | Sin  wt| 

and  the  magnitude  of  the  Fourier  components  of  Sin|wt|? 
is  given  by 


8(2TT)* 

1 -4k2 


if  n is  odd 


= 0 if  n is  even 

where  Tis  the  gamma  function. 

Thus  for  very  large cC only  the  odd  harmonics  are 
present  while  the  even  harmonics  become  insignificant. 


3.2.4  Harmonics  of  the  Flow  Through  a Varying  Orifice 
with  Sinusoidal  Pressure  Difference. 

Another  case  of  interest  is  that  of  an  orifice 
whose  opening  varies  sinusoidally  over  one-half  cycle 
and  is  zero  for  the  other  half  cycle. 


r 
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If  the  pressure  difference  between  the 
orifice  is  sinusoidal  or  predominantly 
pressure  is  in  phase  with  the  opening, 
solution  for  the  harmonic  analysis  can 

follows : 


two  sides  of  the 
sinusoidal  and  the 
a closed  form 
be  obtained  as 


Let  P = P sin  wt 

and  X = X sin  wt  if  sin  wt  >0 


x = o if  sin  wt  -0 


Where  X is  the  linear  orifice  open'  % the 
flow  Q = CpA  ir  or 


q = CpWX  sin  wt 
Where  W is  the  width 


(P  sin  wt)* 
of  the  orifice 


opening  and 


,n  _ a - * (Sin  wt)3/2  if  si»  wt 

A A * 

C WX  P * 

^ s 0 if  sin  wt  iO 


The 


Fourier  series  representation  of  Q is 


Qn  « 4r 


Fa. 


oc 


T + 

n=1 


(A  cos  nwt  + Bn  sin  nwt) 


where 


A = 0 if  n is  odd,  and 


n 
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and  Bn  = 0 if  n is  even,  and 


B = 
n 5 


(n  + 1) 

_ JELtzil 5 L 


V^(i  + t,  j - f) 

where  is  the  Beta  function 


Q = 


.874 

TT 


+ 


cos  wt 

7T 


- .528  sin  2 wt  - .111  cos  3_wt 
-047  sin  4 wt  -.026  cos  5 wt 


The  above  analysis  indicates  that  second  harmonic  is 
significant  (more  than  50$  of  the  first)  while  higher 
harmonics  are  substantially  less  in  magnitude. 


5.3  PARAMETRIC  PERFORMANCE , HARMONIC  ANALYSIS  AND 

TRADE-OFF  DATA  FOR  HYDROACOUSTIC  POWER  OSCILLATORS 


One  form  of  the  hydroacoustic  power  oscillator 
is  shown  in  Figure  3.6  which  utilizes  a variable 
orifice  valve  to  modulate  the  flow  of  sea  water 
from  the  sea  pressure  head  to  an  evacuated  reservoir. 
The  spool  valve  is  driven  with  a nearly  sinusoidal 
force  by  the  oscillator  sections  at  either  end  of 
the  spool. 

The  basic  differential  equations  for  the  power 
oscillators  have  been  programmed  into  a previously 
developed  computer  program  which  solves  sets  of 
coupled,  first-order,  ordinary,  differential  equations. 
This  model  is  then  used  to  simulate  the  harmonic 
content  of  the  pressure  applied  to  the  acoustic 
transmission  line  and  thence  into  the  water.  It  is 
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Hydroacoustic  Power  Oscillator  which  Modulates  the 
if  Sea  Water. 


also  used  to  generate  parametric  data  relating  Input 
and  output  power,  efficiency,  average  flow,  peak 
valve  spool  displacement  and  diameter,  and  pressure 
modulation  coefficient  to  the  acoustic  load  impedance, 
(both  magnitude  and  angle),  fundamental  frequency 
of  oscillation,  and  pressure  head.  This  data  is  used 
to  select  the  optimum  hydroacoustic  power  oscillator 
parameters  for  a wide  range  of  applications.  The 
results  have  been  normalized  with  respect  to  the 
acoustic  load  impedance,  operating  frequency,  and 
acoustic  filter  parameters.  Table  3-1  presents 
the  glossary  of  symbols  and  normalizing  factors  used 
in  this  section. 


3.3.I  Basic  Power  Oscillator  Circuits 

Two  basic  types  of  hydroacoustic  power  oscillators 
are  considered  in  this  report.  The  first  (shown  in 
Figure  3.7a)  is  similar  to  the  device  shown  in  Figure 
3.6  and  is  designated  as  a single-ended  oscillator. 

The  second  (shown  in  Figure  3.8a)  is  designated  as  a 
double-ended  oscillator.  It  will  be  shown  in  this 
section  that  the  performance  of  the  double-ended 
oscillator  can  be  directly  related  to  the  performance 
of  the  single-ended  oscillator. 

Figure  3.7b  shows  an  equivalent  representation 
of  the  single-ended  oscillator  section  of  the  hydro- 
acoustic transducer  shown  in  Figure  3.6.  The  single 
orifice  is  fed  by  sea  water  pressure,  Pg,  through 
a long  tube  of  equivalent  inertance  L and  a cavity 
of  compliance  C.  The  flow  through  the  oscillator 
is  modulated  by  a valve  spool  whose  displacement 
is  given  by  Xj  sin  wt.  A typical  flow  waveform  for 
this  class  of  oscillator  is  shown  in  Figure  3.7c. 
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TABLE  3.1  SYMBOLS  USED  IN  SECTION  3.3 
Hydraulic  Pressures  (psi) 


Ps  Supply  Pressure 

Py  Average  pressure  difference  across  the  orifice 

Pv  = Pg  for  a single-ended  oscillator 

Pv  = Pg/2  for  a double-ended  oscillator 

PL(t)  The  pressure  across  the  load 

PL  = amplitude  of  the  fundamental  of  P^Ct) 

= average  pressure  across  the  load 

?T  = P for  a single-ended  oscillator 
Jj  s 

PL  = Pg/2  for  a double-ended  oscillator 
Hydraulic  Flows  (cis) 


Qg(t)  Flow  through  feed  inertance  L from  supply 

pressure  to  the  valve 

Qg  = average  flow  from  the  supply  pressure  source 

Qy(t)  Flow  through  the  orifice  or  valve 

§y  = average  valve  flow  (equal  to  <5g) 


QT ( t ) Flow  through  the  load 

Jj 

QL(t)  = PL(t)/RL 

= average  flow  through  load  = 0) 

Ql  = amplitude  of  the  fundamental  of  QL(t) 
Amplitude  of  harmonic  of  a flow 

Valve  Parameters 

Ly  Metering  diameter  of  the  valve 


TABLE  3.1  (continued) 


Porting  fraction 


Porting  length  (K  =fTDyPp) 


Xs(t) 


Valve  displacement  Xg(t)  = XT  sin  wt 
XT  = amplitude  of  valve  displacement 
w = angular  driving  frequency 
f = driving  frequency  (f  = 2Tfu>) 


Discharge  coefficient  of  valve  (C  = 92  for 


sea  water) 


- ciW57p 


Cjj  = contraction  coefficient 


f - density  of  fluid 


Total  flow  inertance  in  series  with  valve. 


Ly  = 2 f/K  + + Lj 


Lyg  = Upstream  flow  inertance 


= Downstream  flow  inertance 


Load  Parameters 


Load  resistance 


Inlet  path  flow  inertance  (L  = pi/A) 

j>  - fluid  density 

fi  = length  of  flow  path 


A = cross-section  area  of  flow  path 


II 
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TABLE  3.1  (continued) 


It  ! 


Load  cavity  compliance  (C  = V/jp  c.  ) 
V = load  cavity  volume 
f = fluid  density 
C = speed  of  sound  in  fluid 


Load  Q factor  (parallel  R-L-C  circuit) 

Q = H 0 = Vote 


<t> 


Phase  angle  of  load  on  valve 


Tuned  frequency  of  load  circuit 


f = 2T TU>0 

0 


U>„  = 1/(10)* 


Performance  Parameters 


W, 


Input  power  (V^  = PgQy/8.86  Watts) 


W. 


F 


Output  power  (WQ  = PjQj/17.72  Watts) 
Efficiency  ( \ = WQ/Wi) 


Normalizing  factor  for  flows 
= Cq^**V  XT 


F, 


W 


RL 


Normalizing  factor  for  power 
FW  = CqtfDV  XT 

Normalizing  factor  for  load  resistance 
FRL  ^ / ^Cq"^DV  XT^ 
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TABLE  3.1  (continued) 


Normalized  Variables 


Normalized  frequency 


Kw  = w/wQ  = f/f( 


Normalized  load  resistance 


KRL  " rl/frl 


Normalized  input  power 


Win  * VFW 


Normalized  average  flow  through  valve 


^Vn  " ^v/fq 


Modulation  Coefficient  (X=  P^/P-^) 


0t  = ?T/P_  for  single-ended  oscillator 

L o 


OC  = 2Pl/Ps  for  double-ended  oscillator 


Harmonic  distortion  index 


EDI  = 


(.4  ■*’)/(&  •»’) 


■ i v r m jr 


E 


K 


I 


RETURN 
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a.  Schematic  Diagram 


El, 


T 


i. 


h1 


-1, 


{ (t 


b.  Equivalent  Representation 


- C 


c.  Typical  Waveform  for  Flow  Through  the  Valve 


Figure  3.7  Basic  single  Ended  Hydroacoustic  Power 
Oscillator 
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The  flow  finally  goes  into  a low  pressure  reservoir 
(assumed  here  to  have  a constant  pressure  of  zero). 

Figure  3.8a  is  a sketch  of  a symmetrical,  double- 
ended,  power  oscillator  with  two  variable  orifices  that 
alternately  modulate  the  flow  into  and  out  of  the  load 
cavity.  Figure  3.8b  is  the  analogous  representation  of 
Figure  3.8a.  In  this  device,  the  flow  is  caused  by  a 
pressure  differential  between  the  supply  pressure  and  the 
average  pressure  at  load  driving  one  half  of  the  cycle. 
During  the  other  half,  the  flow  occurs  from  the  load 
cavity  to  the  return  pressure  line  through  another 
orifice.  Since,  on  a steady  basis,  there  is  no  flow  through 
the  load  and  the  average  flow  into  the  load  cavity  during 
first  half  cycle  must  be  equal  to  average  outflow  during 
the  other  half,  the  mean  differential  pressures  at  load 
relative  to  supply  pressure  and  to  the  return  pressure 
must  ce  equal  as  seen  in  Figure  3.8c,  the  flow  waveform 
for  a double-ended  oscillator  is  similar  to  that  for  the 
single-ended  case  for  the  first  half  cycle  the  flow 
waveform  is  repeated,  but  with  the  opposite  sign. 

The  performance  parameters  for  the  double-ended 
power  oscillator  can  be  directly  related  to  the  para- 
meters for  the  single-ended  case.  When  the  two 
configurations  are  operated  at  the  same  output  power 
level  with  the  same  supply  pressure,  the  average  dif- 
ferential pressure  across  the  valve  for  the  double-ended 
case  is  Pg/2  and  Pg  for  the  single-ended  case.  This 
leads  to  the  following  relationships  between  single 
and  double-ended  parameters. 
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c.  Typisal  Waveform  for  Flow  Through  the  Valve 


Figrue  3.3  Basic  Double  Ended  Hydroacoustic  Power 
Oscillator 


(Dy  Xt)d  = y[2  (Dv  XT)S 

<Vd  = 

($L)D  = 2 (Ql)s 

(Wo^D  = (Wo^5 
( i ) D = ^Vs 

= Ws 

( d = (°Og 

(Rl>d  = <rlV4 

(pS^D  = (PS^S 

Comparing  the  Fourier  series  -expansion  cf  (QyCt))^  to 
that  of  (Qv(t))s*  shows  that  the  even  harmonics  in 
the  double-ended  case  are  zero,  and  the  odd  harmonics 
(including  the  first)  are  twice  those  of  the  single- 
ended  case. 


(Qv(t))D  = f(wt)  wt  = 0 tolT 

=-f  (wt-fT)  wt  = 7T  to  2 TT  . 
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3.3.2  Normalizing  Factors 


The  performance  of  many  different  power  oscillators 
operating  into  many  different  loads  can  he  expressed 
in  several  graphs  using  nondimensional  quantities 
for  the  output  parameters.  In  this  report  the  output 
parameters  of  interest  are  the  average  flow  through  the 
valve,  Sv»  (which  is  equal  to  the  average  flow  from  the 
power  source,  ^g)*  the  input  power,  W^;  the  output 
power,  WQ;  the  efficiency,.^;  (0\ie  equal  to  WQ/Wi); 
the  amplitude  of  the  fundamental  harmonic  of  the  pressure 
at  the  load,  PL;  and  the  harmonic  content  of  the  pres- 
sure at  the  load.  The  normalizing  factors  for  the  first 
three  parameters  (Qy  VL  , and  WQ)  are  given  in  Table  3«1 
along  with  the  normalizing  factors  for  the  load  and  oper- 
ating frequency.  Efficiency  is  already  normalized  so 
no  factor  is  needed  in  this  case.  The  first  harmonic  of 
load  pressure  is  normalized  by  dividing  it  by  the  average 
pressure  at  the  load  and  calling  the  result  the  modulation 
coefficient CL. 


The  harmonic  content  is  normally  measured  by 
a factor  called  Total  Harmonic  Distortion: 


The  typical  harmonic  analyzer  which  is  used  to 
estimate  this  distortion  normally  measures  another 
factor  defined  below. 


Harmonic  distortion  = 
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This  factor  is  close  to  total  harmonic  distortion 
if  the  second  and  higher  harmonics  are  not  large 
compared  to  the  first  harmonic*  In  this  report, 
the  square  of  the  above  quantity  is  used  which  is  called 
here  the  harmonic  distortion  index  (HDI)  i.e; 

t 4 

HDI  = 


k = 1 


The  first  nine  harmonics  are  included  in  HDI. 


If  the  driving  frequency  is  different  from  the 
resonant  frequency  of  the  filter  circuit,  the  load 
will  not  be  truly  resistive  but  will  have  a reactive 
component  whose  sign  depends  on  the  frequency.  It 
can  be  easily  shown  that  for  a parallel  F.-L-C  circuit, 


the  phase  angle  # is  given  by 

i - (f/f0)2 


tan0=  TfTj; 


The  results  from  using  three  phase  angles  (0°,  45‘ 
and  50°)  are  presented  in  Section  3*3.4. 


3 Mathematical  Model 


The  behavior  of  the  single-ended  hydroacoustic 
power  oscillator  shown  in  Figure  3.7  can  be  described 
by  the  following  differential  equations. 


dQV  = 


"V 


PQv 


L + 2KL,  W ~ 


2*  Vxs 


For  Xs>  0 
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dQv 

cTT’ 


= 0 


‘v  = o 


dQS 

~ns 


dPL 

"It 


_1_ 

L 


For  X $0 
S 


P - P 
S Ii 


= r [v 


«v  - Vel 


■] 


The  variables  Qy,  Qg  and  P^  are  all  zero  at 

tine  equal  to  zero.  These  equations  are  solved 
using  a previously  developed  computer  program  which 
performs  a numerical  integration  of  coupled,  non- 
linear, first  order  differential  equations.  The 
time  integration  is  continued  until  the  transient 
part  of  the  solution  has  died  out  and  only  the  steady- 
state  solution  remains. 


3.3.4.  Presentation  of  Results 


Figures  3.9  a,  b,and  c are  plots  of  the  modulation 
factory  and  the  harmonic  distortion  index  KDI  as  a 
function  of  the  normalized  resistance  K^t  for  different 
values  of  Q (1.0,  3.0  and  5.0  respectively).  These 
plots  have  three  curves  each  for  the  two  variables  ot 
and  HDI,  representing  three  typical  values  of  phase 
angle  0of  the  load.  The  modulation  factor  * in  every 
case  goes  up  sharply  until  reaches  about  8.0  and 
then  tends  to  level  off.  Also  there  is  no  remarkable 
variation  in «< from  one  case  to  another.  The 
magnitude  of  the  first  harmonic  of  the  pressure  in 
the  load,  P^,  is  given  byOCPg.  The  presence  of 
higher  harmonics  is  measured  by  tha  harmonic  distortion 
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index  as  previously  defined.  Significant  variation 
in  this  parameter  is  expected  as  Q,  a measure  of 
filtering  characteristics,  is  increased.  Figure  3.9a 
indicates  that  harmonic  distortion  is  greatly  increased 
as  the  load  phase  angle  increases  when  the  driving 
frequency  shifts  away  from  the  resonance  frequency 
of  the  load  circuit.  A similar  trend  also  is 
indicated  by  Figures  3.9  b and  c although  the  relative 
magnitude  of  the  difference  becomes  less  predominant 
as  Kbt  increases.  The  value  of  HDI  decreases  sharply 

ftll 

as  Q is  raised  from  1 to  3.  A minor  reduction  results 
as  Q increases  from  3 to  5.  Also,  for  every  case, 
the  value  of  HDI  goes  up  slowly  up  to  values  of  KRL 
from  6.0  to  8.0,  then  shows  an  almost  abrupt  rise  in 
rate  of  increase  and  then  tends  to  level  off  for 
values  of  KR^  around  18.  There  is  also  a slight 
reduction  (more  noticeable  when  Q is  equal  to  1 ) 
in  HDI  as  KR^  varies  from  1.0  to  5.0. 

Figures  3.10  a,  b,  and  c are  plots  of  efficiency 
^ as  a function  of  normalized  resistance  K^. 
Efficiency,  lik efl(,  increases  rapidly  for  values  of 

between  0.0  and  6.0.  For  Q equal  to  1.0,  the 
efficiency  increases  monotonically  for  increasing 
Kr^  and  shows  a considerable  decline  as0increases 
from  0°  to  45°  and  then  to  60°.  For  Q equal  to  3.0 
or  5.0,  the  same  trend  exists  for  values  of 
the  efficiency  becomes  greater  for  higher  values  of  <t> 
as  indicated  by  Figures  3.10  b and  c. 

Figures  3.11  a,  b,  and  c form  a set  of  plots  of 
Win  and  WQn  as  a function  of  KRL  where  Win  and  WQn 
are  the  normalized  input  and  output  power  respectively. 
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Normalized  Values  of  Input  Power,  Output  Power  and  Average  Flow 
as  a Function  of  K^yfor  Different  Values  of0.  Q=5«0 


Since  input  power  is  equal  to  the  average  flow  through 
the  valve  times  the  constant  supply  pressure,  Pgj  the 
input  power  curves,  when  properly  scaled  (^yn=Win/17.72) 
also  represent  the  normalized  average  flow  through 

the  valve.  The  input  power  decreases  monotonically  as 
Kj^  increases  and  is  generally  higher  for  larger  The 
output  power  peaks  at  KRL  approxim  tely  equal  to  3.0  or 
4.0  and  then  gradually  decreases  for  higher  values  of  KDT 

rvJj  • 

The  maximum  output  power  for  a given  Q is  largest  for  0 
= 0°  and  decreases  as ^increases. 


Figures  3.12  a,  b,  and  c are  mostly  the  replots 
of  previous  curves  to  assess  the  impact  of  increasing 
Q for  a given  value  of  0 on  the  various  output 
parameters  likefy*,  and  HDI.  Figure  3.12a 
indicates  that  the  value  of  0 has  little  effect  on  ol. 
up  to  KR^  equal  to  6.0  for  Q equal  to  1 up  to  10.0 
whereas  there  is  slight  variation  for  higher  values 
of  KRI*  The  maximum  divergence  (on  the  order  of  15$) 
occurs  at  the  highest  value  of  Kj^  shown.  The  efficiency 
also  shows  almost  no  variation  up  to  Kj^  equal  to  5.0 
whereas  beyond  that,  the  curves  diverge  rapidly  and 
at  Krl  equal  to  20.0,  n varies  from  .455  for  Q equal 
to  5.0  to  0.72  for  Q equal  to  1.0.  As  expected  the  most 
significant  variation  in  HDI  occurs  as  Q is  increased. 

For  all  values  of  Q the  HDI  curves  dip  down  slightly 
as  ^R1  ^ncreases  from  1.0  to  about  4.0  and  then  rise  for 
increasing  KRL*  Thus  it  would  seem  that  optimum  value 
°f  Krl  for  reducing  harmonic  distortion  for  a‘  given  Q is 
about  4.0. 

Figures  3.13  a,  b,  and  c are  the  plots  of  HDI 
for  the  flow  through  the  valve  (as  compared  to  previous 
plots  of  HDI  for  flow  through  the  load,  i.e.,  flow 
after  filtering)  for  Q equal  to  1.0,  3.0,  and  5.0 
with  varying  <p.  These  curves  are  drawn  to  examine 


Figure  3.12a  Variation  of  as  a Function  of  KRL  for  Different  Values  of  Q 


Variation  of  /K  as  a Function  of  KRL  for  Different  Values  of  Q 


Figure  3.12c  Variation  of  HDI  an  a Function  of  KpT  for  Different  Value 
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Figure  3»13c  Variation  of  HDI  of  the  Flow  through  the  Valve  as  a Function  of 
Kdt  for  Different  Values  of  Q.  0 = 60° 


the  harmonic  distortion  before  filtering  and  to 
measure  the  effectiveness  of  the  designed  filter. 

Figures  3.14  a,  b,  and  c are  the  plots  of  HDI 
as  a function  of  KRL  for  the  double-ended  device. 
These  values  are  computed  from  the  data  of  single 
orifice  device  as  explained  before.  The  harmonic 
distortion  in  general  is  correspondingly  less  than 
that  for  the  single  orifice  device)  the  amount  of 
reduction  depends  upon  other  parameters  like  Q 
and  Krl. 

3.3.5  Trade-Off  Considerations 

The  design  objectives  of  a hydroacoustic 
transducer-amplifier  are  to  achieve  high  efficiency, 
high  output  power  and  low  harmonic  distortion. 
Fortunately,  much  trade-off  is  not  required  as 
optimum  conditions  exist  for  a certain  set  of  para- 
meters. After  a review  of  the  description  of  the 
curves  in  the  previous  section,  the  following 
points  are  worth  noting: 

1)  The  optimum  value  of  Kj^  lies  between  3.0  and 
5.0  with  slight  variation  depending  on  other 
parameters.  In  this  range,  the  efficiency  is 
highest  or  close  to  it,  the.  harmonic  distortion 
is  minimum  and  the  output  power  is  at  its  peak. 
The  device  should  be  operated  at  its  resonant 
frequency.  The  value  of  Q does  not  make  sign- 
ificant effect  onilor\in  this  range  but  affects 
HDI  considerably.  A higher  value  of  Q would 


T?ifmre  3 14a  Variation  of  HDI  for  the  Double 
Figure  3.14a  v r Different  Values  of  Q. 
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be  desirable  and  an  optimum  value  will  be 
dictated  by  any  physical  constraints  (size  of 
cavities,  etc.)  and  the  permissible  harmonic 
distortion.  As  seen  in  Figure  3.12c,  the  HDI 
is  reduced  dramatically  when  Q varies  from  1.0 
to  3.C  and  only  slightly  when  Q is  raised  from 
3.0  to  5.0,  The  reduction  i3  still  smaller  as 
Q varies  from  5.0  to  10.0.'  Thus,  if  the  device 
is  „o  operate  at  a constant  or  almost  constant 
driving  frequency,  an  optimum  or  near  optimum 
design  can  be  easily  achieved. 

2)  One  trade-off  would  be  an  apparent  choice  between 
the  single  and  double-ended  device.  The  double- 
ended  device  leads  to  a low  value  of  HDI; 
comparison  of  HDI  curves  in  Figure  3.12c  (single- 
ended)  and  in  Figure  3.14a  (double-ended)  in- 
dicates roughly  50$  reduction  in  the  value  of  HDI. 
At  Kri  about  4.0  and 0 equal  to  0°  this  reduction 
in  an  already  small  index  is  not  very  significant. 
However,  at  other  values  of  #and  KRL,  the 
reduced  index  may  require  attention  as  a design 
factor.  It  should  also  be  noted  that  for  a 

gi  .’en  power  rating,  the  double-ended  device 
requires  a larger  stroke  (by  a factor  of ^2) 
and  two  metering  edges. 

3)  If  the  driving  frequency  f varies  over  a certain 
range,  the  phase  angle  0 of  the  load  will  not  be 
zero  and  will  vary  according  to  the  value  of  f 
and  Q.  If  Q is  large,  a given  change  in  frequency 
will  cause  a greater  change  in  phase  angle 
relative  to  the  change  if  Q was  small  for  example; 
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4569 


Q = 1 § ip  - 60° , w/wQ  = . 

Q = 5 <p>=  60° , w/wQ  = .8417 

As  noted  before,  a higher  value  of^will  lead 
to  greater  harmonic  distortion  in  general.  Thus, 
depending  upon  the  value  of  f a desirable  value 
of  Q should  be  selected  from  the  curves.  Such 
a value  should  be  corresponding  to  highland  OL 
a small  value  of  HDI  over  the  range  of  frequency- 
change  . 

3.3.6  Other  Areas  of  Investigation 


Exploratory  investigation  was  made  in  other 
fields  related  to  acoustic  flows  as  originally  proposed 
while  stress  was  placed  on  the  areas  which  were  of 
direct  and  immediate  benefit. 


High  amplitude  acoustic  flows  in  pipes  and 
through  orifices  are  known  to  be  turbulent  although 
a criterion  for  determining  the  existence  of  tur- 
bulence (like  critical  Reynolds  number  for  steady 
flows  through  pipes)  is  not  available  for  acoustic 
flows.  In  the  absence  of  a detailed  mathematical 
analysis,  the  Reynolds  number  for  such  flows  may 
be  defined  by  Rg  = Ud  where  U is  the  average  root 
mean  square  velocit?  U = ^ /*  „ (y)y  dy>  „here 


T 


R is  the  radius  of  the  pipe,  d is  the  diameter, 
u (r,t)  defines  the  velocity  field  in  the  pipe, 
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r represents  the  radial  axis,  T is  the  averaging  time 
period,  (much  longer  than  the  life  time  of  large 
eddies),  and  O is  the  kinematic  viscosity. 

For  plane  waves  of  the  type 

u = u exp  (-CU  X exp  (i£Kr  X-ct^) 

Where  «;is  the  attenuation  factor,  <*r  is  the  wave 

, A 

number,  Uis  the  amplitude  of  the  plane  wave  at  X equal 
to  zerOj and  c is  the  phase  velocity. 

A simpler  way  to  define  Reynolds  number  for  acoustic 

flows  is  Re  = 

V 

Neglecting  the  contribution  by  heat  conduction 
Kirchoff's  formula  for  damping  factor 


.M-is  the  coefficient  of  viscosity.  Since  turbulence 
increases  the  dissipation  losses,  its  effect  can 
be  accounted  for  by  replacing  the  coefficient  of 
viscosity^* by  the  eddy  viscosity  coefficient  £ . 

Lacking  the  mathematical  tools  to  handle  the  theory 
to  turbulence,  the  use  of  eddy  viscosity  has  been 
a familiar  way  to  estimate  the  additional  losses. 

However,  there  is  no  way  to  estimate  analytically. 

The  value  of  £ is  generally  determined  from  experi- 
mental data  regarding  steady  flow  in  pipes  etc. 

Binder  (*)  used  a similar  technique  for  measuring 
the  damping  of  large  amplitude  vibrations  of  a fluid 
in  a pipe,  and  by  adjusting  the  value  of £ found  a 
fair  agreement  with  experiment.  A more  sophisticated 

♦Binder,  R.C.,  "The  Damping  of  Large  Amplitude  Vibrations 
of  a Fluid  in  a Pipe",  Journal  of  the  Acoustic  Society 
of  America;  P41 , Volume  15,  1943. 
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handle  on  the  turbulence  losses  will  require  an 


intensive  look  into  experiments  on  purely  acoustic 
and  mixed  (steady  and  acoustic)  flows  whex’e  loss 
measurements  have  been  made.  4n  approximate  analy- 
tical Eodel  is  also  probably  feasible  on  lines  of 
the  turbulent  flow  through  jets  and  pipes  where 
eddy  viscosity  is  related  to  other  parameters  like 
velocity.  In  the  mean  time,  a crude  determination 


may  be  made  by  the  approach  defined  above. 
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